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The purpose of this paper is to study the optimal performance for an irreversible quantum Stirling cooler with 
heat leak and other irreversible losses. The relationship between the optimal cooling load and the coefficient of 
performance (COP) for the quantum Stirling cooler is derived. The maximum cooling load and the corresponding 
COP, as well as the maximum coefficient of performance and the corresponding cooling load are obtained. The 
experimental observation for the optimal region is provided. 
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Introduction 

The finite time thermodynamics, founded by 
Novikov 111 , Curzon and Ahlbom 121 , is a powerful tool for 
the investigation of real thermal devices. It is a discipline 
to study the optimum effect of thermal processes 
originating in finite time and the optimum performance 
of thermal devices and systems operating in finite time. 
In recent years, finite time thermodynamics has made 
considerable progress 13 ' 41 . Geva and Kosloff 15 81 and Wu, 
Chen and Sun 14,9 121 introduced a dynamical semi-group 
approach of quantum mechanics and non-equilibrium 
statistical theory into the finite time thermodynamics. It 
has a wide range of application for performance analysis 
of practical engineering cycles. 

Stirling (ST) coolers have the advantages of high 
efficiency, low vibration levels, low noise, long service 
life, simple structure, etc. Therefore, they have been 
under investigation for the last several years for use in 
infrared devices and in cryogenics. Several authors have 
studied the finite time thermodynamics performance of 
the ST cooler. Wu, Chen and Sun 1 " 1 have discussed the 
optimal performance for forward and reverse 
endoreversible quantum Stirling (QST) cycle. Based on 
the above work, the optimum performance for an 


irreversible QST cooler with heat leak and other 
irreversible losses is undertaken by using the 
non-equilibrium statistical mechanics and the finite time 
thermodynamics. The fundamental relationship between 
the cooling load R and the coefficient of performance 
(COP) £ is derived in this paper. Planck’s constant 
/i = 6.63xlO -34 (J-s) and h =h/2 n = l.05xl0- 34 (J s) 
in the International System, but h -1 is set in the 
calculation of the paper (It is equivalent to use different 
unit system). 

Irreversible QST Cooler 

The QST cooler cycle is made of two isothermal 
branches connected by two regenerating branches. The 
following assumptions are made for the QST cooler 
model in this paper: 

(1) The working fluid of the cycle consists of many 
non-interacting harmonic oscillators at the atomic level. 
The population of the oscillator n at thermal equilibrium 
can be obtained from the Bose-Einstein distribution 1 ’ 31 


exp(/3co) -1 

where (0 >0 is the frequency of the oscillator, (5 =1/T , 
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and T is the internal temperature of the working fluid. 

(2) The cooler operates between an isothermal heat 
sink B h at temperature T h and an isothermal heat source 
B c at temperature T c . The heat source and sink are 
infinitely large and their internal relaxation is very strong. 
Thus the disturbance of the working fluid on the source 
or the sink can be ignored. Therefore the source B c and 
the sink B h are assumed to be in thermal equilibrium. 
There is heat leak Q r between the sink B h and the source 
B c . It is the result of the coupling action between B h and 
B c by the medium of the cooler. 

(3) The internal temperatures of the working fluid 
are T\ and T 2 in the two isothermal processes. We assume 
T 2 >T e (T e is the critical temperature) such that the 
Bose-Einstein condensation of the working fluid will not 
be considered. The second law of thermodynamics 
requires /? 2 > p c > P h> P i, where /?,- =1/7) (i = h, c, 
1 and 2). 

(4) In the two regenerating processes, the 
frequencies of the oscillator are constant®. The 
frequencies of the oscillator in the two regenerating 
processes connecting the two isothermal processes are 
0) i and CO 2 with (0 i > CO 2 . The QST cycle plotted in 
the (n, CO ) plane is show in Fig.l. Assuming that the 
internal temperature of the working fluid in the 
regenerating processes changes uniformly with respect to 
time 1 " 1 , we have 

(2 ) 

dt dt k 

where k = k\ k 2 , k\ is the time coefficient of the 
regeneration which is independent of temperature but 
dependent on the property of regenerative materials, and 
k 2 is the delay-time coefficient which is dependent on the 
phase relation between the piston of the compressor and 
the free displacer, inside which the regenerator is. The 
signs “+” and correspond to the processes 4—1 and 
2—3, respectively. 

(5) The irreversible cold losses inside the cooler 
make real cold capacity far less than that obtained by an 
endoreversible theoretical analysis. The coefficient of 
cold losses (f) i is defined as 

I Qi^iQi (3) 

Where Q 2 is the cold capacity obtained by the 
endoreversible theory 141 , and 2 Q, is the sum of the cold 
losses inside the cooler besides the leak Q r . 

Cycle Period 

The Hamiltonian of the harmonic oscillator system S, 
H s , is described in the following form 151 

H s = oW = taa+a (4) 


where N is the number operator, a ~ and a are the 
Bosonic creation and annihilation operators { N = a + a, 
[a,a+]=a a + -a + d= 1}. 

In the two isothermal processes, the working fluid 
system S becomes an open system owing to the coupling 
action between S and the heat sink B h (or the source B c ). 
Using the generalized master equation of an open system, 
and in the Heisenberg picture, one obtains the motion 
equation of the operator x in the system S |4- ' J ' 111 

dx 

- = i(L s +T rB L sB p B )x- 

dt 

I{ {(Oji) + [x,Qj]Qi - (coy )- Qi[x,Qj]} (5) 
ij 

where L s is the Liouville’s operator of the system S, 
T rB represents the trace of the operator of the heat 
reservoir ( B h or /?,.), L sB represents the coupling action 
between the working fluid S and the heat reservoir (B h or 
B c ), f> B is the density matrix of the reservoir (B h or B c ), 
(i (Oji) + and ( coy)~ are constant related to the matrix 
elements of the density matrix p B , <2; and Qj are 
the operators of the working fluid S and [,] is the 
quantum Posson’s bracket. 

Solving equation (5) by setting x = N, Q\ =d+, 
Q 2 -a and neglecting the disturbance of the working 
fluid on the reservoir (L s +T rB L sB p B = 0) gives 

. _dn_d<N>_ 
dt dt 

a{ [exp()3 jco) - l]n -1} exp (f/3 jco) (6) 

where a and f are both constants, j = h and c correspond 
to the processes 1—2 and 3—4, respectively. 

Substituting equation (1) into equation (6), 
expanding the time integral to the second order in ($j(0 
at the limit (n » 1 and hfijCO «1) gives 



where t\ and f 2 are the times of the isothermal processes 
1—2 and 3—4, A = (1/g>! -l/a 2 )la > 0 is a constant 
which is independent of temperature but dependent on 
the quantum performance of the working fluid, pjco «1 
in the paper is corresponding to hfijco «1 in the 
international system. 

From equation (2), the times h and t 4 of the 
regeneration processes 2—3 and 4 — 1 are 


( 9 ) 
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The cycle period time, t = t l +t 2 +t 3 +t 4 is 


T = A 


1 

Pk-Pl 


1 

+ Pi-Pc 




( 10 ) 


where v =2k/A. 

Cooling Load and COP 

According to the quantum mechanics theory tl3! , the 
energy of the working fluid is 

U =< H s >-< fil/V > (11) 

From equation (11), the rate of change of energy is 
obtained by 


and creation operators, and the subscript h,c stand for the 
sink and the source, respectively. Setting n c as the 
thermal phonon number of the source gives 

h c = -c{ [exp (P h co c ) -1 ]n c - 1}exp (gp h co c ) (20) 

where c and g are constants. 

Combining equation (20) and n =-—-—r— 

exp03 c <o c )-l 

gives 

Qr = ®cK = 


, [ exp^mj-l 
c Lexp08 c ft) c )-l 



exp (gp h <»c) 


( 21 ) 


U =(b<N > +co < N >= (bn + con (12) 

Comparing equation (12) with time derivative of the first 
law of thermodynamics, we obtain the instantaneous heat 
flow q and the power P in the following equations 

q =(On (13) 

P = -dm (14) 

Substituting equation (1) into equation (13) and 
expanding the heat quantity integral to first order in 
PjO) (i = h, 1, c and 2) at the limit (n » 1 and /3,cu « 
1) gives 


2] = P todn = 

X 

(15) 

Jn i 

a-A 

A 

Pi 

(16) 


where Q { and Qi are the heat exchange quantities in 
the processes 1—2 and 3—4, the constant 
X = lnfni, / <u 2 ) > 0 is the natural logarithm of the 
frequency ratio. 

The real cooling quantity Q' 2 in a cycle is 

Q'i= Qi-IcQi (17) 

P 2 

where <f> 2 =1-0, can be obtained by the method of 
fitting the experimental data. 

The Hamiltonians of the heat sink B h and the heat 
source B c are respectively 

Hh - bh (18) 

H c =(O c b c + b c (19) 

where co k , co c are the frequencies of the thermal phonon 
of the sink and the source, and b, b + are the annihilation 


where Q r = — is the rate of heat leak from the sink 
T 

to the source. 

Because the source and the sink are both infinitely 
large and their internal relaxation is very strong, fi h , 
P c , (o c are all constants. Thus the heat flow Q r is 
constant. 

Combining equations (10), (15), (17) and (21), the 
cooling load R and the COP can be obtained in the 
following equations: 

R _ Q'l-Qr _ 

T 


<M r yP i 

a P 


+ 


ypi 

yp i -Pc 


+P(;y-l )] -1 -Q r 


( 22 ) 


Qi-Qr _ 0 2 -yP,( 2 r t/x 
(Ql-Qr)-(Qi-Qr) y~<t> 2 


(23) 


where y = P 2 ! P\- 

It may be proved from equations (22) and (23) that 
the R and the COP are of the same extreme point. 
Combining equations (22) and (23) with condition 
( dR /9)8 ]) y =0 (or 9 el 9)3, = 0) yields 


R _ x<t> 2 Mi+py- 

A L y — D 2 


+ P(y-i)]-‘-Gr 


(24) 


e = 


y-0 2 X y-D 2 


+n(y-m 


(25) 


where D = (P C / P k ) 05 • equations (24) and (25) give 
the optimal relation between R and e with parameter 
y. When the heat leak is ignored, i.e., Q r = 0 , the 
optimal relation between the R and the COP can be 
obtained in the following equation: 



196 


Journal of Thermal Science, Vol.ll, No.3,2002 


(26) 

A l + g-i-D 2 £ 

Using equation (24) and the condition dR/dy = 0 
yields expression for the maximum cooling load R m , 


ft ___ Q 

" A{[l + (l + /t°-5)Z>] 2 -//} ^ 


The COP corresponding to R m is 

£ <t>2^-AQ r [(\ + D + Dfi 0 - 5 ) 2 -^] 

0 A[D(\ + D)/n-°- 5 + D 1 -<f> 2 ] 

When Q r - 0, equations (27) and (28) become 

D __^02_ 

Am 

A[l + Z) + D/x 0 - 5 ) 2 -n] 


cn----- JU 

D(l + D)/x -°- 5 + D 2 ~<t > 2 

Using equation (25) and the condition de / 3y = 0, 
the maximum COP ( E m ) and the corresponding cooling 
load R 0 can be expressed as 

= —^—(02 ~^f~Wy 0 -!)+— 1 + 7 ^]} ( 31 ) 
y 0 -<t>2 X yo-D 2 

Ro = + + Myo - D]- 1 - 6, (32) 


A y 0 -£> 2 


where 


yo =[«2 +(^2 2 -ai22 3 ) a5 ]/ai 



a! = AQrP-fa IX-AQ r ^l X-<p 2 + AQ r {\ + D) 1 1X Fig-2 R!R m versus e characteristics: the effects of n 


a 2 = D 2 {AQ r n^ 2 lX-AQ r ^lX-<f> 2 ) 
a 2 =D 2 a 2 -AQ r <p 2 D 2 (\ + D) 2 IX 

Discussion 

The optimum cooling load versus the COP 
characteristics depend on X / A , tj > 2 , Q r , /z and )3 C / f5 h . 
The general relationship of the optimal cooling load 
versus the COP for a QST cooler is shown in Fig.2 and 
Fig.3. The curve is a loop-shaped one, thus we obtain the 
optimal region of a QST cooler, i.e. 

R 0 <R<R m ,e 0 <£<£ m (33) 



where R 0 , R m and e 0 , e m are determined by 
equation (32), (27) and (28), (31). 


Fig.3 RIR m versus e characteristics: the effects of 4> 2 
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Conclusion 

The finite time thermodynamic model of a QST 
cooler has been established in this paper. It provides the 
application foundation of the non-equilibrium statistical 
mechanics to the practical ST cycles. The solution of the 
generalized quantum master equation of a thermal 
system is obtained in the Heisenberg picture. It makes 
the quantum picture of the coupling action between the 
working fluid and the source (or sink) more clear and 
simple. 

The optimization region (or criteria) for QST cooler 
is obtained in the paper. The relation between the cooling 
load R and the COP, the maximum cooling load R m 
and the corresponding COP, as well as the maximum 
COP and the corresponding cooling load for an 
irreversible QST cooler with heat leak and other 
irreversible losses are derived. The analysis provides a 
new theoretical basis for the performance evaluation and 
improvement of real ST cooler. 
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